r •• 

7 


N AS A T EC H N I CAL 


REPORT 


' Os ' ; ‘ -V h / 

. co , ,\- 

.S"'/ eo' : ' >■ ’ ’>■ 



m NASA TR R-339 


Cs L 


‘r' ' 4 s \ 


■}y 




)-K. V, v i v *> 


'si* v,^(; 




;i - . V. 


'■Vi- ‘ ,rA - ^ 

V-'' w ') . 


j-jt f,AS\ 

V- ' ‘ - ; '• 


_ v- r,' ' v 
- O- ,,-V 
a «... . , ;• 


.r’« \ ,V 7 ^' ^ • V?^' "'f 

.A..-r=> ,, ' ,; •• , 0 - 1 .. . \ 


) ; . 




\ ' ^ -■*' -V 


OF ELECTROMAGNETIC WAVES 
S OBLIQUELY INCIDENT UPON 
A PLASMA SLAB ’ h ' C fir- 


C" • s » r C' ■ 


;V -'ir- T ' 
x ’ o-iu'fx ; . < 




■; _ V ^ : Sc C; J I 

LQA2I COPY.: raTlIRN. . 

- Am iWWJL); “■> ; 

' XIRJLHTO APB-, H MEX » 

?. r V -a ■■ ... i- ,' r •, 
'.;• •-; f, J “-K+ ( S 

■ ' ' V - •• S .,- 


-n / •• •. 


>.V'- ,'W4f 

. ;■ ■, r; ■ r , 

. ■' ' i 

-. ► > w., < - 


... • • v . : - > ' . 

t -‘j y:^.J'c;. vV' 

‘ . / •’ s 4 ' . • 


byf Calvin .CfT ' Swift: Vc \ 
Langley Research Center 

''\: : O Sr/-''. < • \ •’* ^ - \ " i ■ i*."i 


0A; 


-l ■= * 






' T ' r 1 \ L ' 


. V. J 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION WASHINGTON, D. C, ;• APRIL 1970 


)V.„ •. v 


'U •• 4 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



QDbfl3Dfl 


1. Report No. 2. Government Accession No. 

NASA TR R-339 

3. Recipient's Catalog No. 

4. Title and Subtitle 

A THEORETICAL INVESTIGATION OF ELECTROMAGNETIC 
WAVES OBLIQUELY INCIDENT UPON A PLASMA SLAB 

5. Report Date 

April 1970 

6. Performing Organization Code 

7. Author(s) 

Calvin T. Swift 

8. Performing Organization Report No. 

L-6938 

10. Work Unit No. 

129-01-22-05-23 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23365 

11. Contract or Grant No. 

*13. Type of Report and Period Covered 

Technical Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 

The information presented herein was offered as a thesis in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy in Physics, College of William and 
Mary, Williamsburg, Virginia, August 1969. 


16. Abstract 

The problem of an electromagnetic wave obliquely incident upon a plasma slab is con- 
sidered as a boundary-value problem by use of a self-consistent solution of the coupled 
linearized Vlasov and Maxwell equations. Power reflection, transmission, and absorption 
coefficients are derived under the assumption that all particles undergo specular reflection 
at the surfaces of the plasma slab. Although the analysis is valid for arbitrary slab thickness, 
computational results are presented for slabs which are thin when compared with a wavelength. 
The results show that a series of resonances occur which are attributed to the finite tempera- 
ture of the plasma. The results further show that the resonances are Landau damped as the 
thermal velocity of the plasma electrons increases. It is shown that similar resonances can 
be predicted from the coupled linearized hydrodynamic Maxwell equations; however, as is well 
known, such a model does not predict Landau damping. The effects of a finite collision fre- 
quency are then included by means of a simple Bhatnagar-Gross-Krook (BGK) collision term. 
The numerical computations vividly indicate that the resonances undergo severe damping for 
extremely small ratios of the collision frequency to the signal frequency. 

Finally, the plasma capacitor problem is considered, and the results indicate that the 
longitudinal resonances have characteristics very similar to those of the plane-wave 
resonances. 


17. Key Words (Suggested by Author(s) J 

Electromagnetic waves incident upon plasma 
slab, Kinetic plasma effects, Landau damping, 
Compressible plasma, Longitudinal plasma j 
oscillations 


18. Distribution Statement 

Unclassified - Unlimited 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 


21. No. of Pages 

56 


22. Price* 

$3.00 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 





CONTENTS 


SUMMARY 1 

INTRODUCTION 1 

SYMBOLS 4 

INTERACTION OF A PLANE WAVE WITH A UNIFORM PLASMA SLAB 9 

Interaction of a Plane Wave With a Uniform Cold Plasma Slab 10 

Fluid Description of Plane Wave Problem 13 

Direct Solution of Plane Wave Problem Using Linearized Vlasov Equation .... 16 

Kinetic Results for Electric Field Perpendicular to Plane of Incidence 24 

Relationship to Plasma Capacitor Problem 26 

Numerical Results 28 

RESULTS AND DISCUSSION 37 

CONCLUDING REMARKS 39 

APPENDIX A - FURTHER PROPERTIES OF THE DISPERSION INTEGRALS .... 41 

APPENDIX B - INTERACTION OF A PLANE WAVE 

WITH A PLASMA HALF SPACE 47 

APPENDIX C - SOME COMMENTS ON THE COMPUTATIONAL PROCEDURE ... 51 

REFERENCES 53 


iii 



A THEORETICAL INVESTIGATION OF ELECTROMAGNETIC WAVES 
OBLIQUELY INCIDENT UPON A PLASMA SLAB* 

By Calvin T. Swift 
Langley Research Center 

SUMMARY 

The problem of an electromagnetic wave obliquely incident upon a plasma slab is 
considered as a boundary-value problem by use of a self-consistent solution of the 
coupled linearized Vlasov and Maxwell equations. Power reflection, transmission, and 
absorption coefficients are derived under the assumption that all particles undergo specu- 
lar reflection at the surfaces of the plasma slab. Although the analysis is valid for arbi- 
trary slab thickness, computational results are presented for slabs which are thin when 
compared with a wavelength. The results show that a series of resonances occur which 
are attributed to the finite temperature of the plasma. The results further show that the 
resonances are Landau damped as the thermal velocity of the plasma electrons increases. 
It is shown that similar resonances can be predicted from the coupled linearized hydro- 
dynamic Maxwell equations; however, as is well known, such a model does not predict 
Landau damping. The effects of a finite collision frequency are then included by means 
of a simple Bhatnagar -Gross -Krook (BGK) collision term. The numerical computations 
vividly indicate that the resonances undergo severe damping for extremely small ratios 
of the collision frequency to the signal frequency. 

Finally, the plasma capacitor problem is considered, and the results indicate that 
the longitudinal resonances have characteristics very similar to those of the plane-wave 
resonances. 


INTRODUCTION 

The interaction of electromagnetic waves with plasmas has been of continuing inter- 
est to those engaged in the study of the ionosphere, of the radar return from meteor trails, 
and of reentry plasma sheaths. In the last category, most of the research emphasis has 
been directed toward the solution of boundary-value problems that vary in complexity 


*The information presented herein was offered as a thesis in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy in Physics, College of William 
and Mary, Williamsburg, Virginia, August 1969. 



from relatively simple plane wave interactions (ref. 1) to rather complicated ones 
involving antennas under "cold" plasmas (ref. 2) and compressible plasmas (ref. 3). 

A cold plasma is defined here as one in which the electron thermal velocity is zero, 
and thus the plasma behaves as an incompressible fluid which exerts no pressure. One of 
the primary shortcomings of this model is that no mechanism is provided in which to 
excite a spectrum of longitudinal plasma waves. The importance of these longitudinal 
oscillations lies in the fact that they have been observed experimentally in the laboratory, 
as far back as 1931 (ref. 4) and later in connection with the radar scattering from cylin- 
drical plasma columns (ref. 5). It was observed that the radar return consisted of a 
series of resonances, the characteristics of which are shown in figure 1. The interesting 

feature of these resonances is that cold 
plasma theory predicts only the main 
resonance which, for the cylindrical 
column, occurs at o> p = ^ 20 ) where cu 
is the frequency of propagation and u> p 
is the plasma frequency. (For a plane 
slab, this resonance occurs when w = u> p .^ 
Much more comprehensive and convincing 
experiments were initiated in 1957 by 
Dattner (ref. 6). His experiment consisted 
of placing a small cylindrical discharge 
tube into the wall of a waveguide with the 
electric field perpendicular to the tube and 
monitoring the reflection coefficient as a 
function of increasing discharge tube cur- 
rent. Because of the thoroughness of the experiment, there is no doubt that the secondary 
resonances exist, and they have since that time been termed Tonks -Dattner resonances. 
Although excellent experimental results were available, a satisfactory analytical explana- 
tion for these resonances did not appear until the classic work of Parker, Nickel, and 
Gould (ref. 7) was published. Upon applying a fluid model of the plasma, they were able 
to conclude that a spectrum of resonances were generated at the frequencies 

9 O 3KTp 
P m 

where T e is the electron temperature, K is the Boltzmann constant, m is the elec- 
tron mass, and k is a wave number which depends on the radius of the plasma column. 
These frequencies correspond to longitudinal plasma oscillations, which couple more 



Figure I.- Radar return from a cylindrical plasma column. 
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strongly to the electromagnetic wave because the phase velocity of such waves is of the 
order of the thermal velocity of the electrons in the plasma. 

One of the shortcomings of an approach based on the linearized fluid equations lies 
in the fact that for finite thermal velocity, the fluid equations are valid only for 10 ~ cop 
(ref. 7). Analytical results based on such a model should therefore become less inaccu- 
rate as the ratio of Vp/co decreases. Another shortcoming is the failure of the fluid 
equations to predict Landau damping (ref. 8). The proper description requires a detailed 
solution of the Vlasov equation, where the Vlasov equation is identified as the collisionless 
Boltzmann equation to be solved self-consistently with Maxwell's equations. The impor- 
tance of the Landau damping lies in the fact that the collisionless damping should be pro- 
nounced as the ratio of the thermal velocity to the phase velocity of the longitudinal wave 
increases. It ydll be shown that this ratio increases as the order of the resonance 
increases, and could account for the damping of the secondary resonances shown in fig- 
ure 1. In order to determine how the widths of the resonances at half -maximum behave 
in detail as a function of all the parameters, it is necessary to solve a boundary -value 
problem. The model considered herein consists of a plane wave obliquely incident upon 
a plasma slab with the electric vector polarized in the plane of incidence so that longi- 
tudinal plasma oscillations are excited. This particular problem also has applications 
to the study of antennas under reentry plasmas, because the radiation characteristics of 
such antennas can be described by a spectrum of plane waves. (See ref. 9.) The kinetic 
treatment of this problem has previously been considered by Hinton (ref. 10) and by 
Bowman and Weston (ref. 11) in the United States and by Kondratenko and Miroshnichenko 
(ref. 12) in the Soviet Union. Hinton (ref. 10) solved the problem by expressing the cur- 
rents as integrals over particle orbits. This procedure is equivalent to solving the 
Vlasov equation. The approach, however, requires several ponderous perturbation 
expansions and leads to an integral equation solution of the problem. Bowman and Weston 
(ref. 11), on the other hand, used the singular eigenfunction techniques of Case (ref. 13), 
Shure (ref. 14), Felderhof (ref. 15), and Van Kampen (ref. 16) to obtain solutions to the 
Vlasov -Maxwell equations. The disadvantage of this approach is that analytical and 
numerical results appear to be rather difficult to obtain. Kondratenko and Miroshnichenko 
(ref. 12) published an excellent and concise piece of work. Proceeding as Landau (ref. 8) 
did for the half-space problem, they used an integrating factor to solve the Vlasov equa- 
tion. This procedure resulted in a solution in the form of an integral equation which was 
reduced by means of a Fourier series. The treatment presented herein differs from 
theirs, largely in the initial formulation procedure, although the mathematical results are 
the same. 

In none of these papers were numerical results presented. In fact, the only compu- 
tations which have appeared were done by Melnyk (ref. 17), who considered a plasma the 
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equilibrium statistics of which are governed by degenerate Fermi -Dirac statistics. 
Maxwell -Boltzmann statistics will be considered and the problem will be approached by 
initially assuming specular reflection of electrons at the plasma boundaries. This pro- 
cedure, as will be shown, automatically allows an immediate choice of a Fourier series 
representation of the problem. This procedure does not lead to a solution expressed in 
the form of an integral equation. It is in this way that the formulation herein differs from 
that in reference 12. The usual electromagnetic boundary conditions are used in connec- 
tion with the boundary condition of specular reflection. The reflection, transmission, 
and absorption coefficients are then solved for and calculated as functions of the plasma 
electron density and thermal velocity for a slab which is thin compared with a free-space 
wavelength and for zero collision frequency. A series of resonances, that is, peaks in the 
reflection coefficient, occur which exhibit features of the Tonks-Dattner resonances, and 
which become Landau damped as the thermal velocity of the plasma increases. The 
reflection coefficient described by a continuous fluid model of a plasma is also computed; 
similar resonances are noted except that they are not Landau damped. 

A kinetic analysis of the plasma capacitor (ref. 18) is included to strengthen the 
physical deduction concerning the predominance of longitudinal oscillations in the plane - 
wave solution. The results show that the plasma capacitor, which contains only a longi- 
tudinal electric field, resonates at precisely the same slab thickness, plasma frequency, 
thermal velocity, and propagating frequency as those for the plane wave interacting with 
the slab. These resonances are more conventionally defined in the sense that a peak in 
resistance and a zero in reactance are noted at the resonant frequency. 

Finally, a finite collision frequency is considered by using a simple Bhatnagar- 
Gross-Krook (BGK) collision term (ref. 19), and for purposes of nomenclature, the kinetic 
equation will be referred to as the Vlasov equation. The results show that the higher 
order resonances are completely damped out at such a small value of the ratio of colli- 
sion to propagating frequency that laboratory reproduction of such resonances would be 
difficult to achieve at normal radio and microwave frequencies. It is concluded that 
although the present model exhibits some characteristics of the Tonks-Dattner reso- 
nances on a qualitative basis, the detailed structure of the resonances is influenced by 
another mechanism, probably the inhomogeneity of the plasma. 

SYMBOLS 


A power absorption coefficient 


a 


speed of sound in electron gas, 



4 



B magnetic flux density 

b = (a> - koV z sin 6)/v x 
C capacitance 

c speed of light 

D dielectric displacement 

E electric field intensity 

E 0 amplitude of incident electric electric field intensity 

e electronic charge 

F 0 normalized probability function (assumed to be Gaussian in this paper) 

F+ = f + + f“ 

F" = f + - f“ 

f + distribution function for particles having velocity v x greater than zero 

f” distribution function for particles having velocity v x less than zero 

f distribution function (no velocity half-plane restrictions) 

Gj,G 2 impedance coefficients for parallel polarization defined in equation (15) 

G]^ G 2 ^ impedance coefficients for perpendicular polarization defined in equation (79) 
H magnetic field intensity 

H 0 amplitude of incident magnetic field intensity 

Hj,H 2 * Fourier transforms of two variables in velocity space, defined in 
equations (62) 
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I/A 


current per unit area 


i-fi 

^ 1^11)^21^31’^ 41’^ 51 dispersion integrals defined by equation (Al) 

j current density 

K Boltzmann constant 

k wave number 

k 0 wave number in free space 

kp x-component of wave number in a cold plasma; kp = k 0 

Wpy'w < cos 6 and k p = ik Q \jsin^6- ~ for Wp/cd > cos 6 

k u wave number of longitudinal plasma oscillations; k u = ICq 

for — ^ € -~ < sin 6 
a/c 

L thickness of plasma slab 

l Fourier indices (l = 0,1,2, 3,. . .) 

m mass of electron 

n electron number density 

P pressure 

R reflection coefficient or resistance, as appropriate 

T transmission coefficient 

T e electron temperature 
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— > sin 6 and k u = ik^/sin^S 
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Mill ■IIMMUMII 


II I I mill 


immi mu 



time 


t 

u 

u X> u yj u Z 

V 

V 

v ph 

v-p 

X 

X 0 

x,y,z 

x' 

z 

y 

bj 

e 

e o 

? 

<r 

e 


fluid velocity 

unit vectors in three principal directions 
voltage 

particle velocity 
phase velocity, u>/k 
root-mean-square particle velocity, 
reactance 

reactance of air-filled parallel -plate capacitor, l/o>C 

Cartesian coordinates 

dummy variable for integration 

impedance 

constant in adiabatic equation, Pn - ^ = Constant 
Kronecker delta, 0 when i ^ j and 1 when i = j 
permittivity 

permittivity of free space 

argument of dispersion function, = r I 

U /, k 0 2 L 2 sln 2 e 

r * l V 

value of £ for l = 0 

value of £ for case of a nonzero collision frequency 
angle of incidence 



K 


wave number, In/L 


A 

A + 

** o 
v 

a 

u> 

cop 


V 

[] 


velocity transform variable 


= 1 - 




permeability of free space 
collision frequency 
conductivity 

angular frequency of propagation 
angular plasma frequency 
vector 
tensor 

absolute value 
del operator 
matrix 

Fourier transform 


Subscripts: 

x,y,z denotes vector components along the principal directions 

l denotes Fourier components 

0 denotes first-order or unperturbed quantity, unless otherwise specified 

1 denotes second-order or perturbed quantity, unless otherwise specified 
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denotes component of reflection or transmission coefficient perpendicular to 
the plane of incidence 

denotes component of reflection or transmission coefficient parallel to the 
plane of incidence 

INTERACTION OF A PLANE WAVE WITH A 
UNIFORM PLASMA SLAB 


Figure 2 shows the geometry of the problem. A plane wave is incident upon a 
plasma slab with the electric vector polarized in the plane of incidence. The incident 
electromagnetic wave is assumed to have a harmonic dependence of the form 

E = E 0 exp ifk Q x cos 8 + IJqZ sin 8 - cot) 


The faces of the plasma slab are x = 0 and x = L; the plane of incidence is the 
x,z plane; and the angle of incidence is 8. Here k D is the free-space propagation con- 
stant k Q = o>/c and co is the angular frequency of the incident wave. The case where 

the electric vector of the incident wave is per- 
pendicular to the plane of incidence is discussed 
in an analogous manner in a later section. Kinetic 
effects, however, depend upon the ratio of thermal 
velocity to the phase velocity of the plasma waves 
involved. This ratio is appreciable only for longi- 
tudinal plasma waves. These waves, however, do 
not couple to incident electromagnetic waves 
polarized perpendicular to the plane of incidence. 
The reflection coefficient for the case where the 
electric vector is parallel to the plane of inci- 
dence is discussed, first for the linearized cold 
plasma model, then for the linearized fluid 

model, and finally for the linearized Vlasov equation (with a BGK collision term). The 
equations describing the plasma in each case are solved self-consistently with Maxwell's 
equations. In each case, the tangential field components at the left of the slab (x < 0) are 
given by 



Figure 2.- Geometry of a plane wave obliquely 
incident upon a plasma slab. 


Hy _ H 0 (e ik ° x cos 0 + Re -ikoX cos 0 ) sin 6 “ wt ) 


( 1 ) 
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( 2 ) 


E z = cos 6 - Re -i ^o x c ° s 0)^5 C os 6 e 1 ^ 2 sin 6 ~^ 


where R is the complex reflection coefficient for the magnetic field and H 0 is the 
magnetic field amplitude of the incident wave. MKSA units are used throughout the paper. 

To the right of the slab (x > L), 


H = HoTe^ 031 cos ^e^^o 2 sin 6 -a,t ) 


( 3 ) 


and 

E z = -H 0 T 

Here T is the complex transmission coefficient for the magnetic field. The boundary 
conditions across the surfaces x = 0 and x = L require that the z -dependence of the 
fields within the slab be the same as those outside; therefore, the fields inside the slab 

are of the form [ E,H] = [ E(x),H(x)j e^o 2 s * n ® ~ As such, the exponential depend- 
ence e^o 2 S3n ® ^ need not explicitly appear in any of the subsequent expressions. 


COS £ p i(koZ sin 6 - u>t) 


cos 6 


( 4 ) 


Interaction of a Plane Wave With a Uniform Cold Plasma Slab 

If the plasma is cold, the random velocity of the free electrons is assumed to be 
zero, and the dielectric constant of the plasma can be determined without resorting to 
kinetic theory. The equations of motion of a free electron interacting with an electro- 
magnetic wave are solved in order to deduce the polarization per particle; from the 
polarization the following expression for the relative dielectric constant of the plasma is 
obtained: 




( 5 ) 


2 _ n o e " 


where u>p is the plasma frequency, = — — (n 0 is the electron density, m is the 

** * IUCq ' 

electron mass, e is the electron charge, and e Q is the permittivity of free space) and 

The plasma frequency a>p and the 


v is the collision frequency for momentum transfer, 
collision frequency v are assumed to be constant. 
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The solution to this problem appears in reference 20 but not in a form which will 
be useful when these results are compared with those obtained for the fluid and kinetic 
models. 

In order to develop the desired solution the Maxwell curl equations 

VXH=^ = -io>e E (6) 

ot 


VXE = -^ = i^ 0 H 


(7) 


are used. Since the plasma is nonparamagnetic, the permeability of the plasma /i is 
assumed to be that of free space p. 0 , and the dielectric constant e is given by equa- 
tion (5). By using equations (6) and (7), H = Hy(x) Uy in the slab (uy is the unit vectoT 
in the y-direction) can formally be derived. The result is 


Hy(x) = 


Hy(L) sin kpx - Hy(0) sin kp(x - L) 
sin kpL 


( 8 ) 


where Hy(0) and Hy(L) are the values of Hy at x = 0 and x = L, respectively, and 

kp = ko^|- - sin^0 for o>p/a> < cos 6 and kp = iko \Jsin%9 - for cop/w > cos 6. 

Reflection, transmission, and absorption coefficients can now be determined from the 
boundary conditions, that is, continuity of tangential E and tangential H at x = 0 and 
x = L. Use of equations (1) to (4), (6), (8), and the boundary conditions leads to the fol- 
lowing relationships: 


(1 + R)H 0 = H y (0) (9) 

H 0 Te ik ° L cos 6 = H y (L) (10) 

-^(1 - R)H 0 cos 0 = -io>iU o |Hy(0) GjL - H y (L) G 2 l] (11) 

-^H 0 Te ik o L cos 6 cos 6 = -iu>M o [H y (0) G 2 L - H y (L) Gil] (12) 

It is important to note here that the functions Gj and G 2 are defined separately for 
(a) the cold plasma model, (b) the fluid model, and (c) the Vlasov model. In this way, a 
single algebraic relationship can be used to solve equations (9) to (12) once Gj and G 2 
are given. For the cold plasma, 
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(13) 


Gl = 


[Vi j_ 
W jL 


cot knL 


c oJ 


g 2 = 


G 1 


cos kpL 


rfg] JL 

^°' fej 


esc kpL 


kpL 


If Z is defined as the surface impedance of the plasma at x = 0, 


Z = - 

cos 9 1 


g 2 " 


Gl - 


i cos 6 

koL 


It follows then that 


(14) 


(15) 


R = 


1 - Z 
1 + Z 


(16) 


and 

_ -iko L cos 9 
(1 + R)G 2 e ^ 

o i cos 9 

1 '~kjr 


The absorption coefficient may be defined as: 

A = l - |r| 2 - |t| 2 


(17) 


(18) 


In the absence of collisions (p- = 0), the absorption coefficient of the cold plasma slab is 
zero. 

In the limit as L approaches °°, only forward-traveling waves exist in the 
plasma slab, and the expression for the surface impedance reduces to 


Z 


^ k P 1 1 

L - ” 00 ko JL cos 6 
e o 


(19) 


Near (a>p/cc>)2 = l and (v/<jS) ~ 0, the dielectric constant approaches zero, and 
Gi and G£ become large; thus, the surface impedance can approximately be written as 
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( 20 ) 


Z 



tan kpL ik 0 L 
kpL cos 8 


which for a thin slab (kpL approaching 0) reduces to 


/ \2 

Z ikpL /kp\ i 

e-4) cos #\k 0 / A 
L-0 ' 0/ 


( 21 ) 


It is seen that the reflection coefficient of even a thin slab, for which v - 0 and e 
approaches 0, should be unity at co = u>p because the impedance becomes infinite. Also 
note that the limit given by equation (21) depends upon the order in which the limits are 
taken with respect to e and L. Naturally, for L s 0, Z goes to one, not zero as 
implied by equation (21). The thin slab is investigated in more detail later. When the 
slab is not thin, from equations (13), (14), and (15) it can be seen that the impedance 
approaches zero when kpL = (l + for l = 0,1,2,. . . . These values are the Fabry- 
Perot resonances which are familiar in optics. Numerical results are given in a sub- 
sequent section. 


Fluid Description of Plane Wave Problem 

The linearized Vlasov equation, with a BGK collision term of the form 
-v(i - f 0 ) = (where f 0 is the unperturbed distribution function and fj is the 

perturbation) and with 8/dt = -iu>, may be written as 


- 8f l 
-icdfl + V • 

9x 


eE 

m 



( 22 ) 


If the zeroeth and first moments of equation (22) are calculated with respect to the 
particle velocity, the following expressions are obtained (see ref. 21, for example) for 
conservation of mass and momentum: 


-iumi + n 0 V • u = 0 


(23) 


(-iw + l^)n 0 u = - V • P (24) 

where u is the fluid velocity, P is the pressure tensor, and n Q and n^ are the 
unperturbed electron density and its perturbation, respectively. 
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The pressure tensor corresponds to terms in the next higher moment, which can be 
eliminated by assuming a scalar pressure and using the equation of state 

P = nKT e (25) 

(K is the Boltzmann constant and T e is the electron temperature) in connection with 
an adiabatic equation 


= Constant (26) 

If the pressure term in equation (24) is linearized, and nj in equation (23) is eliminated 
by using Pj = yn^KTg, the following equations are obtained: 

mn 0 (v - iu>)u = r^eE - VPj (27) 

a 2 mn 0 V • u = io)Pj (28) 

where = yv-p 2 = yKT e ^/ m and Vrp is the thermal velocity. Equations (27) and (28) 
are the same ones Wait (ref. 3) used to derive the reflection coefficients for a uniform 
half-space. The procedure herein is similar to his, except that an additional boundary 
at x = L is included. The electric field and the fluid velocity are related by the 
Maxwell equations 


V X E = iu>/i 0 H 

(29) 

V x H = -iu>e 0 E + n Q eu 

(30) 


when the last term in equation (30) is the macroscopic convection current. 

Equations (27) to (30) can be used to develop wave equations for Pj and H. 

Since the wave equation is a second-order differential equation, a total of four unknown 
coefficients must be determined within the slab (two coefficients for Pj and two for h). 
However, equations (27) and (29) can be used to show that the boundary condition of spec- 
ular reflection, that is, u = 0 at x = 0 and x = L, implies that 


d _El 


9z 


icue 0 9Pj 
8x 
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x=0 

x=L 


x=0 

x=L 



The two unknown coefficients for Pj can therefore be expressed in terms of those for 
H in the slab. This procedure leads to the following solutions: 


Hy(x) 


Hy(L) sin kpX - Hy(0) sin kp(x - L) 
sin kpL 


(31) 


where 


up 2 o)ko sin 6 (H y (0) cos k u (x - L) - H y (L) cos k u x] 




mk, 


sin k u L 


(32) 


u 


k - 2= M i -^- +i 4- k o 2sin2e 


(33) 


Note that the expression for Hy(x) is identical to that obtained for the cold plasma. It 
also follows that 


EZ= ^ 


<Hy(0) 


k p cos k p (x - L) 3~kp 2 sin 2 e cos k u (x - L) 


- sin kpL 


(! + i^ k u sin k u L 


Therefore, if 


and 


+ Hy(L) 




kpcoskpx ^-k 0 2 sin2e cos k u x 
^ sin k p L (l + sin k u L 


W P 2 9 

-Zj- sin^ 6 cos k u L 
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k u L 


(34) 


(35) 
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&>/ ^-kpL sin k p L ^1 + i^)^k u L sin k u L 


(36) 
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an expression for the impedance at x = 0 identical with the cold plasma result (eq. (15)) 
is obtained, but with the functions Gj and G 2 given now by equations (35) and (36). 
When k u = °°, that is, when a^ = 0, Gj and G 2 reduce to those for the cold plasma. 
By inspecting equations (35) and (36), it is noted that for v = 0, Fabry-Perot type 
resonances occur when k u L -Its {l = 0,1,2,. . .) in addition to the cold plasma reso- 
nances e/e c = 0 and kpL = (l + 1/2 )is ( l = 0,1,2,. . .). It is interesting to note at this 
point that the phase velocity can be very low for longitudinal waves, and therefore 
resonances can be expected for slabs which are thin compared with a free-space wave- 
length. If the thermal velocity Vij> = a is small compared with the speed of light, 
these resonances occur when 


k u L = mr ~ 


k o L f t ^2 
v T r ■ w 2 
fy c 


(37) 


or at a phase velocity Vp^ = -^-c. The phase velocity of longitudinal waves therefore 
becomes smaller as the slab dimension decreases and as the order of the resonances 
increases. The quantity y is normally assumed to be three for an electron gas. 

Results of calculations for the reflection and absorption coefficients are presented 
in the section "Numerical Results," and comparisons are made with those obtained from 
the other models of the plasma. 


Direct Solution of Plane Wave Problem 
Using Linearized Vlasov Equation 

In the preceding section, the linearized Vlasov equation (eq. (22)) was written with a 
BGK collision term of the form -v(t - f Q ) = and with 9/3t = -iu> as 


- iuf i + 7 ■ If 


m 9y A 


where f 0 is the unperturbed distribution function, fj is the perturbation of the distri- 
bution function, v is the particle velocity, and x is the particle position. Equa- 
tion (22) is to be solved self-consistently with the Maxwell curl equations: 

V x H = -io>e 0 E - e y fjv dv (38) 

V x E = iwju 0 H (39) 
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where the last term in equation (38) is the density of the convection current. In order to 
proceed, the distribution functions for the velocity half-space v x > 0 and the velocity 
half-space v x < 0 are considered separately. If the former is denoted by fi + (- . . v x . . .) 
and the latter denoted by fj“(. . . -v x . . .), fj + and f^ - satisfy the following equations: 


. 3f i . 

+ v)h + v x + ik 0 v z sin 6 f x 


e_ 

m 



+ 



(40) 


(-iu> + v)i± 


9fi 

v x + ik ° Vz sin 6 fl 



(41) 


If the expressions F + = fj + + fj" and F" = fj + - fj” are introduced, the following 
second-order differential equation in x is obtained for F": 


82f (o> + iv - k 0 v z sin 0)^F" 2ie(u> + iv - k 0 v z sin 6) 9f 0 ^ 2e 9 ^o 9 ®z _ A / 

~dj~ + mv^2 ~§T = ° 

If all particles are specularly reflected at x = 0 and x = L, then F" must vanish at 
x = 0 and x = L. This condition can be satisfied identically by a Fourier sine series 
for F" as a function of x: 


F " = 2 F *('0 sin T £ 

1=1 

with 

F," = f- F" sin 4 s dx 
L Jq l 


(43) 


(44) 


Examining equations (42) and (43) indicates that they imply a Fourier sine series expan- 
sion for E x and a cosine series for E z . 

If Ej x , Ej z , and Hjy are the corresponding Fourier coefficients for E x , E z , 
and H y , 
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Ff = 


2ieo)v x F oi 


]/^V H iy-( 1+ ‘sK 


mVrp^ 


(oj - koVzSin 6 + i^) 2 - 


(45) 


where 


f 0 = n o F o = 


n o e 


-v 2 / 2 v t 2 


^2i7V-p 2 ^ 


3/2 


is the equilibrium distribution function and is assumed to be Maxwellian. The Fourier 
components of F + are derived in a similar manner. Since the x and z components 
of the current density are defined as 


ix - - e jf *y I 0 


*2 *y \„ v x F " *x 


and 


] z = -e dv z dv y V Z F + dv 2 


it follows that the Fourier components of j and j are: 


Vrr ,2 In 

V T 


e o c N 


, _iV we oL( T ,, /^o J 2J 

" v 5> Ztt \ J li E ^X " 1/ £_ C 
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2,.w= , ( [R J 
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_ g>p 0*0 L | u 

] * z ” Vt 2 y e 0 c ^ 


(46) 


where 


J U- ( 1 + i s)f I j’i' dv y^I 0 


v x 2f o ^x 


. \2 (lirv- 


(oo - k Q v z sin 6 + iuj - 


(47) 
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V x v z F o 


/ \2 /iirv x Y 

- k 0 v z sin 8 + ii/j - jr-j 


(48) 


. lA r°°r ^ ^ r 00 v z(“> - k 0 v z Sin 0 )f o dv x 

J3i = 2(1 + l5J ) ^ * y * z —— 2 

(w - ko v z sin 6 + iv ) 2 - Y-l - ] 


(49) 


T „ CC rh . ^ r°° *x v z 2 (“ - k oV z sin 0 )f o 

J 4l - 2 ^ Chr y dv z j 7“ \2 

-OO U l . . n . \2 ( J,rv x\ 

- k o v z sin 0 + w') - Y~l/ 


(50) 


A Fourier analysis of the Maxwell curl equations (eqs. (38) and (39)) gives 


lit 


j; E Zz = -ik D sin 6 E lx + iu,M 0 H z y 


(51) 


E Zx = 


ko sin 6 


u>e r 


H Zy + 

‘•y iu>e. 


1 . 


•*£x 


(52) 


H 


iy 


l^ 2 - k Q 2 sin 2 6 - 


Its 


Its 


H y (0) - (-irHy(L) 


L/2 


= -iko sin 9 j lx - l f ] lz 


(53) 


Therefore, if the modal components of the current (eq. (46)) are substituted into 
equations (52) and (53), three equations which contain the unknowns Ej x , E^ z , and H^y 
are obtained. Solution of these equations yields: 


% 


Its 

fey(0) - (-l)*H y (L) 


“l 

L L/2 J 



2 oi«2o _ (l]L 


COp 

_u_L 

Co 2 

,A 


k^L sin 8 


Ti J 2 1) 


(54) 
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E lx = 


In 

Hy(0) - (-l) Z Hy(L) 



u>e 0 L 

L/2 

\f° v t 2 In c j 


k 0 2 - k 0 2 sin 2 0 - 
_ 

) ' k ° 2 0,2 y T zfa - u sm 6 J 2jj 


- (55) 


E* z = 


ICO/Xq 

H y (0) - (-l)*H y (L) 

(l- S m2e-" p2 " jJ 


L/2 

J 

y 1 sin e w 2 v T 2 5l j 

(■-3H 

k Q 2 - k Q 2 sin 2 6 - | 

Cl) -*> 2 -£ T "if«- l -»«i| 


=r (56) 


By using eqviations (33) to (35), 

oo 

E z (0) = ^ E Zz = -ia>M o [H y (0) Gi L - H y (L) G 2 l] 


(57) 


1=0 


where Gj and G 2 are the functions previously introduced to define the surface imped- 
ance for the cold plasma and fluid models. From equations (54), (56), and (57), it follows 
that Gj and G 2 for the Vlasov model are 



f o 07 \ . .T 

where J 5i is defined as Jgj = ^ sin The reflection and transmission 

coefficients can be derived by substituting Gj and G 2 into equations (15), (16), 
and (17). 
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Before proceeding, it is instructive to examine, in more detail, the dispersion 
integrals which appear in equations (58) and (59). Since it has been assumed that the 
unperturbed distribution function is Maxwellian, it follows that 

I co F °( V X’ V y ,Vz ) = F o( V X> V z) 

and the Vy integration trivially disappears. The integration on v x may be extended 
over the entire real axis by noting that the integrand is an even function of v x . Further- 
more, if a partial fraction expansion is used, the integral becomes 

T U • v\ (T ^x dvz F 0 v x fatA 

3 U = (* + l s)j J taT <60) 

u> - k Q v z sin 8 ^ + iv 

with similar reductions for the other dispersion integrals. 

Use can now be made of Fourier transforms in velocity space (ref. 22) to reduce 
them to single integrals. For example, the zeroeth-order dispersion integral Jj can 
be written as a convolution in the form 


J l ~ 


1 + i 


3SS 


dv x dv z F 0 


llTVy- 

-« w - k^v„ sin 8 — - — + iv 


'o z 


fl + i 1L \ 

= — ^ H 2 (a x ,A z ^ dA x dA z 


(61) 


A OO 

U f 

(2,7)2 Jj 


where A x and A z are the transform variables of the velocity components v x and 
v z , respectively, and Hi (A X ,A Z ^ and H 2 *(a x ,A z ) are given by 


Hl( A x.Az) N 


i(A x v x +A z Vz) 




Leo k o L sin 6 ihv 


■°° Vy ‘ " v 


In In 


dv x dv z 


H2*(A X> A Z ) > ^ 


,-l( A xVx+AzVz) e -(vx 2 +v z 2 )/ 2 v t 2 

2nVrj? 


> (62) 


dv x dv z S e 




T 

2 


The first of equations (62) can be evaluated by using contour integration in the complex v x 
plane. The integration over v x gives 


21 



.vL ^ a k o v z L 
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(a x <( 


The integration over v z gives 


H i - -< 2 »> 2 If e 


. Ax«L(l.ng) 




A x > G 


= 0 


(Ax<C 


Equation (61) therefore reduces to 




1 ‘ ' £ 2 7T 2 / 2 


e e 


and the other dispersion integrals become 


. 2 v T 2 (, ko 2 L 2 sln 2 0'; L 


J2X=£vT 4 ^sine f dA x e 
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> (64) 


(65) 
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L x (68) 
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J 4* = -£v T 2 




ko 2 L 2 

l W 



A x 2k o 2lj2v T 2 
l 2 * 2 


sin z 6 


( 69 ) 


As shown in appendix A, these integrals can further be reduced so that they can be 
expressed in terms of the dispersion function of Fried and Conte (ref. 23). 

It should be emphasized again that Gj and G 2 as previously given specify the 
surface impedance of the plasma, and therefore uniquely specify the reflection, trans- 
mission, and absorption coefficients for the plane wave problem. In appendix B, it is 
shown that results reduce to those of the half-space in the limit as L — °°, as they should. 

As the ratio of v T ^v ph (V ph being the phase velocity of the wave) becomes small, 
the imaginary parts of the dispersion integrals J^ n are negligible, and the real parts of 
the dispersion integrals can be expanded in increasing powers of Vrjyv p ^, as shown in 
appendix A. The results, when applied to equations (68) and (69) are 



g 2 




Resonances occur in Gj and G 2 when 


lir = 


koL 


Vs 


1- 


C0 r 




V T 


w 


(70) 


(71) 


(72) 


which is approximately the same as the hydrodynamic results for co p /ci> ~ 1. With k 
defined as lit/ L, it is noted that the resonances defined by equation (72) occur when 
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( 73 ) 


v ph _ gi _ k 0 L 
v T «v T " vt 
c 

Thus the kinetic effect, that is, Landau damping, should become more pronounced as the 
ratio of VpjyV-j- becomes smaller. Therefore, for fixed v-p/c and fixed kpL, the 
Landau damping should become more severe as the order of the resonance increases. 
Numerical results of the reflection and absorption coefficients are discussed subsequently. 

Kinetic Results for Electric Field Perpendicular 
to Plane of Incidence 

It has been indicated that most of the interesting effects associated with a plane 
wave obliquely incident upon a plasma layer occur with longitudinal plasma waves excited 
in the plasma. As such, the case of the electric vector perpendicular to the plane of inci- 
dence = EyUy in fig. 2 ) is only of secondary interest. However, the results are 
included here for completeness. By proceeding as in the section "Direct Solution of 
Plane Wave Problem Using Linearized Vlasov Equation,” the Fourier coefficients of 
F” are found to be 


-2ia)p 2 kp 2 Fo v yV x H^z 


F Z e u) 


(a> - k 0 v z sin 8 + i&») 2 - 


( 74 ) 


It further follows from the Maxwell curl equations that 


V 2 H z + ko 2 H z = 


( 75 ) 


or 


9 2 H z 
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(k 0 2 - kp 2 sin 2 ^H z = -e J dv v. 
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( 76 ) 
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A Fourier expansion solution of equation (76) gives 
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where 
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If and T' L are defined as the reflection and transmission coefficients for the case 
where the electric field is perpendicular to the plane of incidence, then 

E = E (e ik ° x cos 6 + R i e” ik ° x cos ^ sin 0 “ wt ) 


for x < 0 and 


E y = E 0 T e 


l e i(koX cos 8 + kQZ sin 8 - cut) 


for x > L. If the boundary-value problem is completed, it is found that 


1 + R 
1 - R 


1 

-j- = Z = -ik Q L cos 8 


IV* 


c 1 

G i — i — 

y-i 


Gi + - — z— — 

A koL cos 8 


(79) 


where Z is the surface impedance, and the functions G]/ and G 2 ^" are given by 


G * i ■ &) l ! ; 


° » - sin2e - IM.J . ^ 


CUp^ 


(80) 
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(81) 



Equation (79) is slightly different in form from that previously obtained for the case of 
parallel incidence; this difference occurs because R^", in this case, is the ratio of the 
reflected to the incident electric field. 

Relationship to Plasma Capacitor Problem 

It has already been pointed out that although the ratio of v-p/Vp^ is ordinarily 
negligible for transverse waves, parallel-polarized waves might be expected to show 
interesting kinetic effects, as shown. As further support of this assumption, it is useful 
to compare the plane -wave calculations with corresponding computations in which only 
longitudinal plasma oscillations are excited. The plasma capacitor provides this 
information. 

The problem under consideration consists of a plane -parallel plasma-filled capaci- 
tor, whose plates are located at x = 0 and x = L. An electric field, oscillating at an 
angular frequency w, is applied normal to the plates. In this section it is shown that the 
capacitor exhibits resonance behavior at the same values of koL, (a)p/a>) 2 , and Vp/Vp^. 
In order to effect this behavior, the impedance of the capacitor will be determined, as 
done by Hall (ref. 18) and Shure (ref. 24). 

If E z and v are set equal to zero in equations (40) and (41), and a process 
analogous to that which led to the expression for the x-component of the current density 
(eq. (46)) is followed, the Fourier coefficients of the current density in the capacitor are 
found to be 


JZx " 


ictjp 2 aje 0 L 
v 2 llT 

V T 


11 E lx 


(82) 


which is nothing more than equation (46) with H^y = 0. The continuity equation, relating 
charge density and current density, gives 



- iwe 0 E x = Constant 


(83) 


where I/A is the current per unit area on the plate of the capacitor. A Fourier 
expansion of equation (83) gives 
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Since the voltage between the plates is given by V = - f E x dx, the impedance of the 

J 0 

capacitor is given by 


z !/2 f y i - (-d* 

iwcW A / 2 „ 


(85) 




, con^ L 

where C is the capacitance in the absence of a plasma. If A = 1 > equa- 


Vrp 


tion (85) becomes identical to that obtained by Shure (ref. 24). For « 1, equa- 

tion (85) reduces to 


o 00 
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icoCW Zy 


1 - (-D Z 
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(86) 


An inspection of the denominator of the sum in equation (86) shows that resonances 
occur when 


In = 



(87) 


Equation (87) is identical to equation (72), which defined the resonance condition for a 
plane wave incident upon the slab. Therefore similarities between the plane wave and 
capacitor results should be expected. These similarities are exhibited in the form of 
numerical calculations in the next section. 
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Numerical Results 


In order to delineate clearly the differences between using cold, fluid, and kinetic 
models for a plane electromagnetic wave obliquely incident upon a plasma layer, compu- 
tations of the reflection coefficient for each model were made as a function of ^u>p 
with the following parameters fixed: 


6 = 15° 
k 0 L = 0.1 



and v T /v ph = 0.157. The results are shown in figure 3. Note that for the cold plasma 
(fig. 3(a)), only one resonance occurs, and that is located at o> = Wp. In the fluid limit 

L I ^2 \ 

(fig. 3(b)), a series of resonances occur, corresponding to the zeros of sin\^yl — 2 /’ 

Also note that the main resonance is displaced from ( Wp/co ) 2 = 1 and that each higher 
order resonance is narrower than its predecessor. The results obtained from the Vlasov 
equation are similar to those obtained with the hydrodynamic equations, except that the 


xhe departure becomes more 


resonances do not occur at the same values of 
pronounced with increasing order, but this difference is to be expected since the fluid 
approximation becomes less valid. (For example, compare eqs. (37) and (72).) Also 
note that a resonance is associated with each odd term of the Fourier expansion. 


When the phase velocity of the wave becomes of the same order as the thermal 
velocity of the plasma, electrostatic energy in the wave is absorbed by the electrons and 
converted to kinetic energy. This condition occurs because of the inertial "drag" of the 
electrons, and is the well-known phenomenon of Landau damping. In order to demon- 
strate the effects of Landau damping, the 1=5 resonance was investigated for addi- 
tional values of the parameters. The reflection and absorption coefficients are shown in 
figure 4 for various values of v-p/Vp^, where the abscissa has been grossly exaggerated. 
For v^yvph = 0.157, little Landau damping occurs, as evidenced by a narrow resonance 
having a peak reflection coefficient near unity. However, absorption becomes very 
pronounced with only a small increase in v-p/Vp^. Note that the resonance is heavily 
Landau damped at a ratio of v-p/Vph of less than 0.2. It is somewhat surprising 
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(a) Cold plasma theory. 

(b) Linearized fluid model. 

(c) Model based on the linearized Vlasov equation. 

Figure 3.- Reflection coefficient of a plasma wave incident upon a plasma slab. 
k 0 L = 0.1; V T /C = 10-3; v T /Vp h = 0.157. 



Figure 4.- Reflection and absorption coefficients as a 
function of the plasma electron density for various 
electron thermal velocities. 1=5; k 0 L = 0.1. 


to see so much damping at such a low value of v-p/Vp^, but this condition occurs because 
the width of the resonance is so small. The pertinence of the line width is discussed 
later when collisional damping is considered. The resonant peak also occurs at smaller 
values of (cup/c o^“ as vp /Vp^ increases, as is evidenced in figure 5, where the 
abscissa is again exaggerated. The reflection and absorption coefficients are plotted 
as functions of (u>p for various values of the angle of incidence, in figure 6 with the 
following parameters fixed: 
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Figure 6.- Reflection and absorption coefficients as a function of 
electron density for various angles of incidence. I = 5; k 0 L = 0.1; 
vj/c = 1.16 x 10~3; v-r/Vph = 0.182. 


As the angle of incidence increases from 6 = 5° to 8 = 15°, the peak value of the 
absorption coefficient increases, in large part, because the longitudinal component of the 
electric field increases in proportion to sin 8. As 8 further increases, the transverse 
electromagnetic waves become evanescent within the slab. Since the transverse waves 
and longitudinal waves are coupled, this condition leads to a damping of the resonance. 

I ^p^ o 

±iknx ±ikoxJl ^--sin 8 

Since electromagnetic waves in a plasma propagate as e = e V w } the 

1 

waves become evanescent when 8 = cos _i — , which is about 19° for the case considered 
here. Furthermore, note in figure 6 that there is very little shift in the position of the 
resonance with increasing values of the angle of incidence. 

The reflection and absorption coefficients plotted as a function of /ojJ^ and 
v/u> are shown in figure 7 for the model based on the Vlasov equation, with the following 
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parameters fixed: 


e = 15° 


koL = 0.1 

II . 10-3 

c 

and v T /v ph = 0.157. 
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Figure 7.- Reflection and absorption coefficients as a function of 
electron density for various values of collision frequency for a 
plasma model based on the linearized Vlasov equation with a 
Bhatnagar-Gross-Krook collision term. I = 5; k 0 L = 0.1; 

Vy'c = 10* J ; Vf/Vph = 0.157. 
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For values of v/u less than about 10“®, the collisions do not appreciably influ- 
ence the reflection and absorption coefficients; however, as v/oo increases a couple of 
orders of magnitude, the damping becomes pronounced. Although it may, at first glance, 
seem surprising that such a large effect occurs for v/u as low as 10"^, from figure 7 
it can be seen that the line width is of the order of 10“^. From our general knowledge of 
resonance phenomena, damping is expected to be appreciable whenever the collision fre- 
quency is of the order of the line width. This condition also accounts for the degree of 
Landau damping observed in figure 4. Similar conclusions can be drawn by inspecting the 
fluid results shown in figure 8. 



Figure 8.- Reflection and absorption coefficient as a function of 
electron density with collision frequency as a parameter for a 
linearized fluid equation. I = 5; k 0 L = 0.1; v T /c = I0' 3 ; 
v T /v ph = 0.157. 
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Computer results for the impedance of the plasma capacitor are shown in figure 9, 
when the resistance and reactance (normalized to X Q = 1/wC) are plotted as a function 
of Figure 9(a) gives results for v-j/c = 10“3 (v-p/Vp^ = 0.157); figure 9(b), 

for v-py'c = 1.20 x 10" 3 (vx/Yph = 0.188); and figure 9(c), for v T /c = 1.30 X 10" 3 
(vT^Vph = 0.203). In each case, the normalized slab thickness koL is fixed at 0.1. 
Therefore, except for the angle of incidence, which does not appear in the capacitor 




(a) v T /c = 10' 3 ; v T /v ph = 0.157. 

Figure 9.- The impedance of a plasma capacitor as a function of electron density 
for several values of the thermal velocity. [ = 5; k 0 L = 0.1. 


35 



expression, all pertinent capacitor parameters are the same as those for the l = 5 plane - 
wave resonance. From figure 9 it can be seen that the resonance is very sharp for low 
values of v-p/c, and the width at half-maximum noticeably broadens as the thermal veloc- 
ity increases, as does the reflection coefficient of the plane wave. Note further that if 
figure 9 is compared with figures 4 and 5, the peak of the resonance for the capacitor and 
the obliquely incident parallel-polarized plane wave occurs at precisely the same values 
of the plasma parameters. 



(b) V T /C = 1.20 X 10" 3 ; V T /Vp h = 0.188. 
Figure 9.- Continued. 
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(c) vj/c = 1.30 x 10"3 ; vyvpf-, = 0.2(B. 

Figure 9.- Concluded. 

RESULTS AND DISCUSSION 

The prime objective in the work described herein was to examine analytically and 
numerically the details of the coupling phenomena that occurs between electromagnetic 
waves and longitudinal plasma oscillations. This examination was made in order to 
relate this problem to the Tonks-Dattner resonances, which are known to occur when the 
electromagnetic field is applied to an inhomogeneous plasma in such a way as to couple 
to longitudinal plasma waves. The model of the inhomogeneous plasma used was of the 
simplest type, a thin uniform plasma slab. The specific boundary-value problem con- 
sidered an electromagnetic plane wave obliquely incident upon the slab, and assumed 
specular reflection of plasma electrons from the faces of the slab. By solving for the 
reflection, transmission, and absorption coefficients, the detailed behavior of the reson- 
ances, that is, the shift in resonant frequency and changes in the width at half-maximum 
of the resonances in the reflection coefficient as a function of the plasma parameters 
could be examined. The problem was approached in such a way as to delineate the dif- 
ferences, both from a physical and computational viewpoint, between a cold plasma model. 
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a fluid model, and a model based on the Vlasov equation. The case where the electric 
vector was polarized parallel to the plane of incidence was considered in detail. Only 
this polarization excites longitudinal plasma oscillation, for which the ratio of the thermal 
velocity of the particles to the phase velocity of the wave is not so small as to diminish 
the kinetic effects. The analysis further showed that from the practical viewpoint, only 
the lower order longitudinal plasma resonances (that is, l = 1, 3, and 5) can be supported; 
these resonances are well separated only if the plasma slab is thin compared with the 
wavelength of the incident wave. If the plasma is thick, the resonances become more 
grouped near u> = top, for other parameters fixed, as evidenced from equation (72). The 
cold plasma model predicts for such a slab that only one resonance can be supported, and 
that resonance occurs when the signal frequency equals the plasma frequency. The fluid 
model was found to support a series of resonances; these resonances became narrower as 
the order increased, and are similar in nature to the Tonks-Dattner resonances in the 
sense that the secondary resonances occur at successively lower values of the electron 
density than the main resonance does. Similar resonances were observed when a kinetic 
analysis using the Vlasov equation was undertaken except that electron densities for 
resonance were shifted, and the effects of Landau damping became evident as the thermal 
velocity increased. This result is a manifestation of the fact that the ratio vqy'Vph can 
no longer be considered negligible as the order increases. This damping with increasing 
order is qualitatively consistent with the experimental observations of the behavior of the 
Tonks-Dattner resonances. Such collisionless damping cannot be anticipated from the 
fluid equations. The fluid equations, however, have the advantage of presenting a simple 
physical picture of the standing wave processes that occur within the slab. 

The pertinence of the parallel -plate plasma -filled capacitor problem was demon- 
strated by calculating the impedance of the capacitor as a function of slab thickness, 
plasma frequency, propagating frequency, and thermal velocity. It was found that when 
the slab is driven either by the capacitor or by a plane wave, resonances occur at the 
same values of these parameters, and that Landau damping commences in both problems 
at the same values of the thermal velocity. 

Most of the computational effort pertaining to the Vlasov equation consisted of evalu- 
ating the dispersion function, and using it to determine numerically the series given by 
equations (68) and (69). (See appendix C for details.) This series would, in turn, be used 
to determine reflection, transmission, and absorption coefficients. One of the more sur- 
prising aspects of the computations was that, for the cases considered, the infinite series 
given by equations (68) and (69) can accurately be represented by only two terms in the 
series. These terms consist of the first ( l = 0) term and that odd term ( l = 1,3,5,. . .) which 
corresponds most nearly to resonance for the parameters under consideration. This result 
is to be contrasted with the half-space solution (to which the results reduce as the slab 
dimension goes to infinity) which usually requires more involved computational techniques. 
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The influence of collisional damping of the resonances is then considered. For the 
particular parameters considered, it was observed that unless the ratio of collision 
frequency to signal frequency is less than about 10 _ 4, the third odd resonance is damped 
out. Since a collision frequency ratio of less than 10 "4 for gaseous plasmas is not easy 
to achieve in the laboratory, it is not difficult to conclude that the uniform slab can sup- 
port the experimentally observed Tonks-Dattner resonances. 

It is a temptation to conclude that the inhomogeneity of the plasma does more than 
merely control the spacing of the resonances (ref. 7). The results lead us to believe that 
the inhomogeneity may also tend to broaden the width of the higher order resonances at 
half-maximum, which would make them less sensitive to collisional and Landau damping. 
This broadening could occur either because the inhomogeneity provides a gradual transi- 
tion in the impedance between the plasma and the air interface, or because the inhomo- 
geneity induces a background field which causes resonant trapping of the electrons. An 
investigation of the latter problem (that is, consideration of the background field) is a 
very formidable task, indeed. The description of even the simplest problems has devel- 
oped into enormous and frustrating computer programing projects (refs. 25 and 26), and 
illuminating results are difficult to achieve unless many simplifying assumptions are 
made (ref. 27). However, such results should present some new and rather interesting 
kinetic effects, and the problem is therefore worthwhile to pursue. If the background 
field can be ignored (ref. 28), the problem becomes more tractable (but still numerically 
much more difficult than that for the uniform plasma). Such a solution may be valuable in 
order to determine whether a gradual reduction in electron density at the boundaries will 
broaden the resonances. 

The techniques described here can easily be extended to other problems such as: 

(1) the study of the impedance characteristics of antennas under plasmas (ref. 9), (2) the 
effects of nonspecular reflecting boundaries, and (3) electromagnetic waves obliquely inci- 
dent upon a plasma slab in which a static magnetic field is applied normal to the bound- 
aries. All these problems are of interest and may be approached by extending the tech- 
niques of the current work. 


CONCLUDING REMARKS 

The classical problem of an electromagnetic wave obliquely incident upon a plasma 
slab was investigated by using the coupled Maxwell -Vlasov equations. For the param- 
eters considered, a plot of reflection coefficient against electron density reveals a 
series of resonances which become narrower and Landau damped as the order of the 
resonance increases. This behavior is in qualitative agreement with experimentally 
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observed Tonks-Dattner resonances. However, on a quantitative basis, the computed 
resonances are much closer together than experiment indicates; this effect is probably 
due to the inhomogeneity of the plasma. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 6, 1970. 
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APPENDIX A 


FURTHER PROPERTIES OF THE DISPERSION INTEGRALS 

If the collision frequency v is zero, the various dispersion integrals may be 
defined as follows: 

V 

J 21 
J 31 
J 41 
J 5£ 

where the integration over Vy has been performed to give 

r°° 

J oo dVy F o( Vx ’ V y’ Vz ) = F o( V X ? V z) 
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(Al) 


which is the two-dimensional Maxwellian distribution function. From the integral 
relationships, 


oo oo 

v x F o <&x dv Z = Jj* v z F 0 dv x dv z = 0 
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(A2) 


and 


OO 
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(A3) 


It is possible to show that by multiplying the numerators and denominators of the inte- 
grands of equations (A2) and (A3) by co - kjVj. sin 6 - 


Ittv x 




(A4) 
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J 3Z = sin 6 + ^ J 21 

m 

1 = o>J z - n + k Q sin 0 J 3 ^) 

(A6) 

koL 

J 3f = -Jf sin 6 J ll 

(A7) 


As shown in the text, it is possible to use Fourier transforms in velocity space to 
re-express equation (Al) in terms of single integrals over the transform variable A x . 
The results are: 
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(A8) 

In the limit as v-p approaches 0, the exponential involving v-j- 2 may be expanded 
as a Taylor series to give approximations for the real parts of the dispersion integrals as 
the thermal velocity of the plasma approaches zero. The results are 
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These integrals are now related to the tabulated plasma dispersion function, as 
given by Fried and Conte (ref. 23). Fried and Conte define the dispersion function as 


z©-M 
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which may be recognized as the Hilbert transform of the Gaussian. If a change of vari- 

v„ 

ables is used so that x = — — , with 
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and £ is defined as 
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the velocity transforms may be used to show that 


Z(S) = -pi(?) 
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Further manipulations may be performed to show that 

Ju(C) = - 0;(-r) ^ 2 [ 1 + ? z ^>] ( A17 > 

J 2 lid = 2^ sin e ^l^? 3 [ z (S) - 2? 2 Z(?) - 2?] (A18) 

J 3Z (?) = -2 sin 0 [l + ?Z(0] (A19) 


The identities (A4) and (A5) in connection with equations (A17) to (A19) may be used 
to express J 4 £(£) and «J 5 ^(^) in terms of the dispersion function. For 1 = 0, the 
dispersion integrals become 
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J 10 = J 20 = 0 
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where 


? 0 = — p-^ (A24) 

k 0 v2vj sin 6 

Finally, if H^y had been eliminated in equation (43), the Fourier components of the 
current density could have been expressed in the following form: 
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If the modes of the conductivity tensor a values are defined through Ohm's law, 

= °xxZ%x + °xzl E lz\ 


> (A27) 

Hz ~ °z xl E lx + a zzl E lz 


These expressions agree with those of Hinton, who solved this problem by integrating 
over particle trajectories. 

For a finite collision frequency ratio v/w, the pertinent dispersion integrals are 
given by equations (56), (57), (59), and (60), which are also expressible in terms of the 
dispersion function as 
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where, for a finite collision frequency, 



APPENDIX B 


INTERACTION OF A PLANE WAVE WITH A PLASMA HALF SPACE t 

For the case where the slab thickness L approaches infinity, the problem can be 
solved in the original way that Landau and others have done. The Vlasov equation is 
written in the form: 


3fi p 9f 0 

-iu>fi + v x — + ikofiv z sin 9 - — E s 


3v„ 



(Bl) 


The Vlasov equation is again Fourier analyzed in configuration space, except that 
the half-space requires a superposition of continuous modes rather than discrete modes. 
In other words, fj is expanded as a Fourier integral instead of a Fourier series: 
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A straightforward Fourier analysis of equation (Bl), in connection with the convolution 
theorem gives: 
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If the field components of equation (B3) are expressed in terms of their Fourier 
transforms, 
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for all v x where a bar denotes the Fourier transform. From the definition of the cur- 
rent components, 
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where the J terms in equations (B6) and (B7) are identical to the dispersion integrals 
given in appendix A, Itt/L, being replaced by the continuous wave number k x . The 
currents explicitly appear in the wave equation as follows: 
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The form of equation (B3) suggests the symmetry properties Hy(x) = -Hy(-x), 
j x (x) = -j x (-x), and j z (x) = j z (-x). Applying these properties to the Fourier transform 
(eq. (B8)) gives 
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where a bar over a field component denotes a Fourier transform. The Fourier trans- 
forms of the Maxwell curl equations are: 
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Equations (B9), (BIO), and (Bll) are all that are needed to solve for all the field compo- 
nents. After several ponderous algebraic manipulations, the following results are 
obtained: 
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But, it also follows from equations (54) to (56) that Hj y = -H_£ y , E^ x = -E_^ x and 
Elz - *-lz so that the summation index l can be extended from -°° to +°° to yield 
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and similar expressions for E x and E z . In the limit as L approaches the sum- 
mations convert to integrals so that In/L approaches k x and dk x approaches ir/L. 
In this limit, 


H y = H H iy 
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In which case, the sum of the terms defining H y in equation (54) becomes 
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where H y (L) approaches 0 as L approaches °°. Equation (B19) is the Fourier 
transform of equation (B12). Since equations (B13) and (B14) follow in a similar manner, 
it has been proven that the results for the slab reduce to the half-space results as L 
approaches °°, which was to be shown. 
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SOME COMMENTS ON THE COMPUTATIONAL PROCEDURE 

i 

The programs for all numerical computations were written in Fortran IV for use 
on the Control Data 6600 computer facility at the Langley Research Center. The pro- 
grams applicable to the cold plasma and hydrodynamic models involve only the elementary 
transcendental functions of complex arguments and therefore do not require any detailed 
explanation. The programs relating to the Vlasov equation should, however, be discussed 
in a little more detail. The reflection, transmission, and absorption coefficients were 
determined by evaluating the series defining Gj and G2 (eqs. (58) and (59)), which are 
repeated here again for reference: 



Equations (58) and (59) are then used in a trivial manner in connection with equations (15), 
(16), (17), and (18) to solve for the reflection, transmission, and absorption coefficients. 
The dispersion integrals (the J terms) in equations (58) and (59)) were then expressed 
in terms of the dispersion function Z(£), which for £ < 4 was computed from the differ- 
ential equation: 


Z'(?) = -2(1 + ?Z(d] (Cl) 

under the initial condition that Z(0) = For £ > 4, the asymptotic series was used 

(ref. 23). Calculations were performed only for real values of £ 
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For small v/w, the dispersion function was evaluated by using a Taylor series expansion: 
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where, from equation (A32) £ = i^(l + i^. Equation (C2) therefore reduces to 
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For plasma slabs which are thin compared with a wavelength, the numerical results 
indicated that Gi and G2 converged so rapidly that only the 1 = 0 term and the most 
nearly resonant term contribute to the series for Gj and G2. The resonant term is 

co p 2 T 

that term in odd l for which 1 J17 in the denominator of equations (68) and (69) 
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is a minimum. The dominant 1 = 0 term is simply 
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